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1. Introduction 

The experimental program at RHIC has generated much interest in the properties of the quark- 
gluon plasma (QGP) in the theoretical community, not least amongst lattice gauge theorists. Many 
of the dynamical properties of the QGP can only be determined by studying the relevant spectral 
functions p(ffl,p). These are related to the corresponding euclidean-time correlators via an integral 
equation 

f°°d<a , . , . coshffofr- 1/2771 

The survival of bound charmonium states in Nf = 2 QCD well above the transition temperature [jl|], 
the calculation of dilepton and photon production rates [§] and the calculation of QGP transport 
coefficients ||] are all topics which have been studied on the lattice. Determining a spectral 
function from a correlator, evaluated at a finite number of points in a lattice simulation, is an ill- 
posed problem. However, progress can be made by using the maximum entropy method (MEM) []|] 
to determine the most likely form for the spectral function given both the data and a set of prior 
assumptions for the spectral function. This prior information is encoded in a default model which 
is used in the MEM analysis, and one must be careful that the results obtained from this analysis 
are stable under reasonable variations of the default model. 

We are particularly interested in calculating the hydrodynamical contribution in QCD spec- 
tral functions. Transport coefficients can be obtained from the slope of zero-momentum spectral 
functions at ft) = 0, for example, the electrical conductivity is given by 

p em (ft),0) 

a = hm -, (1.2) 

co->0 6 ft) 

where p em is the spectral function for the conserved vector current. Unfortunately, lattice corre- 
lators are extremely insensitive to the details of spectral functions at small frequencies From 
general considerations, one knows that p (ft),p) ~ ft) for ft) <C T, and in the same regime the kernel 



can be written K(r,co) = (2T/co) + ff{(0/T). Substituting these expressions into Eq. (Fl), one 
sees that the bulk of the contribution of the spectral function in the low frequency regime amounts 
to just a constant shift in the euclidean correlator. Therefore, one expects that the low-frequency 
part of a spectral function is the most difficult to compute using MEM. 

As a first step, we have undertaken a study of meson spectral functions above the transition 
temperature. These functions are obtained from the euclidean correlators G(t,x) = (J(z,x)J^ (0,0)) 
where 7(t,x) are quark field bilinears 

J(T,X) = Xff(T,X)r\tf(T,x), r=l,7 5 ,tf. (1.3) 

We also computed these spectral functions for a number of values of spatial momentum. Chang- 
ing the momentum allows us to investigate whether momentum-dependent features can be found, 
which may be helpful in identifying non-trivial hydrodynamic structure. 

This work is particularly timely because there now exist predictions from AdS/CFT for finite- 
momentum spectral functions in ,JV = 4 Super Yang-Mills theory |8|], and a comparison with 
lattice results is of immediate interest. 
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2. Simulation Details 

This study has been performed on sets of 100 quenched gauge configurations generated with 
a Wilson action on two lattices at temperatures both below and above T c . The parameter values 
for these configurations are a subset of the values used in Ref. [Q. The 'cold' configurations were 
generated on a 48 3 x 24 lattice with /3 = 6.5. The temperature on this lattice was approximately 
160 MeV. The second set of configurations were produced on a 64 3 x 24 lattice with j8 = 7.192. 
In this case T ~ 420 MeV, which is about one and a half times the transition temperature. Stag- 
gered fermions were used for the valence quarks, and propagators were generated at three bare 
quark masses, am = 0.01,0.05,0.125. We are interested in comparing our results to hydrodynamic 
predictions and so we want to compute spectral functions for a number of momenta p/T < 1. How- 
ever, for an isotropic lattice with periodic boundary conditions in the spatial directions, the lowest 
non-zero momentum in units of temperature is p/T = 2nN x /N a , where N T is the number of lattice 
sites in the temporal direction and N a is the extent in the spatial directions. On the hot lattice, this 
corresponds to a lowest non-zero momentum p/T ~ 2.4. To circumvent this problem and reach 



sufficiently small momenta, we use twisted boundary conditions for the valence quarks [1C]. We 
computed propagators at four different twist angles and, using these together with Fourier trans- 
forms, we are able to access ~ 20 low-lying momenta. 

2.1 Z 3 symmetry 

It is well known that the pure gauge action has a global Z3 symmetry. It is left unchanged under 
the multiplication of the temporal links by an element of Z3, the centre group of SU(3). Of course, 
this is a symmetry of the quenched theory only and is broken by the fermion determinant appearing 
in the partition function of full QCD. The Polyakov line acts as an order parameter for the centre 
symmetry and the transition to the deconfined phase corresponds to spontaneous breaking of the Z3 
symmetry. The Z3 vacua can be distinguished by the Polyakov loop expectation value, which can 
be real or complex: (P) = l,exp (±/2tt/3). In Ref. JTT| ] it was shown for staggered fermions that 
the nature of the chiral transition in the real vacuum is markedly different from that of the vacua 
with complex Polyakov loops. Above T c in the real phase the chiral condensate drops to zero with 
vanishing quark mass, signalling the restoration of chiral symmetry. This is akin to what one might 
expect in full QCD. In the calculation of the quark propagators we chose our link variables such that 
the phase of the Polyakov loop is real, and observe a similar signal of chiral symmetry restoration: 
above T c for the lightest quark mass m/T = 0.24 the pseudoscalar and scalar correlators are seen 
to coincide. 

2.2 MEM analysis 

For Kogut-Susskind fermions, the lattice correlator receives a contribution from an unwanted 
staggered partner. In this case the expression for the correlator in Eq. < \l.l\ ) becomes 

G(T ' p) = 2 i, ^(^Mp^pM-^p^p)]- (3-D 

To compute p(o,p), we perform independent MEM analyses on odd and even time slices and 
add the results to obtain the desired spectral function. In doing this, one has to be careful that the 
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Figure 1: The plot on the left shows the reconstructed correlator in the vector channel plotted along with 
the original Monte Carlo data points. The right-hand side shows zero-momentum spectral functions in the 
vector channel both above and below T c . 



correlators defined on alternate time slices yield sensible, positive semi-definite spectral functions 
that can be handled by MEM. The main disadvantage of using staggered fermions compared to 
other fermion formalisms is that only half the available time slices are used in each MEM analysis. 
All the results presented here have been obtained using the standard Bryan analysis, although we 
have checked this approach against the classic and historic MEM methods. The default model 
initially used in our analysis is motivated by the expected behaviour of the spectral function in the 
continuum at asymptotically large frequencies Pdefauit {oo) ~ oo 2 . 



3. Results 



As a first check, we confirmed that it was possible to reconstruct the euclidean correlators from 
the results of the MEM analysis. The output for p(ft),p) and p(ft),p) were used to evaluate G(r,p) 
by performing a simple numerical sum which approximated the integral in Eq. (2.1). An example 
of this is shown in the left-hand plot in Figure [j] This plots a reconstructed vector correlator, for the 
lightest quark mass, at vanishing spatial momentum below T c , against the original numerical data 
points. The Monte Carlo results include a sum over the three vector polarisations. The agreement 
between the original data and the reconstructed data points is excellent, and we find this to be true 
in general provided that the data is of high quality with sufficiently small error bars. As one might 
expect, when the data is noisier the match between the reconstructed correlator and the original 
data is less precise. 

In the confined phase below T c , the spectral functions are peaked at energies which correspond 
to hadron resonances. Above T c the spectral functions look very different, particularly in the small 
frequency region. For example, the contribution from Landau damping means that these spectral 
functions are non-zero for all oo > 0. The right-hand side of Figure [l] shows vector meson spectral 
functions at zero spatial momentum above and below T c . These have been rescaled by a factor 
of oo 2 to better reveal their detailed structure. In this plot, the melting of the hadronic state is 
evident. We also note that the spectral functions in both the hot and the cold phase show a large 
bump, centred at co/T ~ 30, which can attributed to lattice artifacts [12, 13 ]. 
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Figure || shows the vector spectral function for varying values of spatial momentum. Below T c 
the position of the spectral function peak shifts to higher energies with increasing spatial momen- 
tum, and the corresponding dispersion relation yields a speed of light which is close to unity. Note 
also that, although we are in the cold phase, the temporal extent of the lattice (~ 1 .2 fm) is not large 
enough to allow us to compute the bound-state energies using conventional lattice techniques, i.e. 
by fitting to the exponential decay of the correlator at sufficiently large times. 

Above T c , in the small co region, there are very clear differences between spectral functions 
with different values of spatial momentum. In Figure |3| there appears to be a threshold energy 
which increases with increasing momentum. The spectral functions are odd functions of co and 
the peak observed at co ~ is an MEM artifact. Determining the precise origin of this artifact 
will require further analysis. For the purpose of this discussion, we adopt a naive approach and 
simply disregard the unphysical peak. The result of this is again shown in Figure [|. The spectral 
density has a bump at small values of co which dies away with increasing momenta. This structure 
resembles the non-trivial structure expected from hydrodynamics and similar behaviour has been 
observed in an AdS/CFT calculation of spectral functions of the stress-energy tensor in JV = 4 
SYM [0]. 
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Figure 2: Vector spectral functions at the lightest quark mass for a range of momenta. 



To test the reliability of these results, we repeated our analysis of the correlators with another 
default model. In this case we chose Pdefauit(w) ~ tanh(ft)/<5) (8 2 + co 2 ) so that for co <C 8 the 
default model varies linearly with co rather than quadratically, and it again scales as co 2 for co » 8. 
For this analysis we took 8 = 0.1. A comparison of the zero-momentum vector spectral functions 
obtained using the original and new default model is shown on the left-hand plot in Figure |]. 
Here, one can see that for co/T > 5 the spectral functions are indistinguishable. Of course, there 
is no reason why they should differ, given that the two default models are virtually identical in 
this region. However in the small co region considerable differences appear. A close-up of this 
is also shown in Figure |[ Once again, we have removed the unphysical spike at the origin. We 
note that at low frequencies the spectral functions obtained using different default models do not 
agree quantitatively, although this difference may be consistent with the uncertainty in the MEM 
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Figure 4: A comparison of zero-momentum vector spectral functions in the hot phase determined using 
different MEM default models as described in the text. The errorbar in the left-most plot gives an indication 
of the uncertainty in the MEM fits. 



fit indicated by the errorbar in the left-hand figure. On the other hand, one might also argue that 
the original default model is completely at odds with the known behaviour of spectral functions for 
small co, and is not an appropriate default model for this regime. We also point out that, because 
of staggering in our correlators, we have used no more than twelve time slices in our analysis 
and it is possible that the results will converge as more time slices are used in the fit. Given the 
precise agreement between the reconstructed correlators and the Monte Carlo data shown earlier, 
this discrepancy highlights the fact that lattice correlators are remarkably insensitive to the details of 
the spectral function at low frequencies. This result indicates how difficult a precise determination 
of transport coefficients will be. 



4. Conclusions/Outlook 



We have computed meson spectral functions at non-zero momentum using the maximum en- 
tropy method. Preliminary results have been obtained both in the cold confined phase and in the hot 
phase above T c . Twisted boundary conditions have been employed to access a range of low-lying 
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momenta in order to look for non-trivial hydrodynamic structure. The maximum entropy method 
can clearly distinguish between spectral functions with differing momenta. Qualitatively, the spec- 
tral functions show little variation under a change of default model. For example in the confined 
phase the positions of peaks which identify hadronic bound-states appear to be stable. However, a 
quantitative determination of the spectral functions at low frequencies is still lacking, and further 
work is required. 

The results presented here were obtained using staggered valence quark propagators. This is 
not the optimal choice for the MEM analysis. We are currently generating quark propagators using 
the non-perturbatively improved clover action [|14|]. This will allow us to include all time slices of 
the resulting correlators in the MEM analysis, which we hope will benefit the analysis. 

Ultimately, it may be that a much more sophisticated approach, including, for example, the 
use of highly anisotropic lattices Jl|,[i5t], is required to precisely determine dynamical properties of 
QCD above T c . However, the exploratory study described here certainly represents a progression 
towards this goal. 
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